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i^ ' Abstract 

00 ; 

^^ ■ The divergences of both the length and time scales, at the 

pg i magnetization- reversal transition in Ising model under a pulsed field, 

^P I have been studied in the linearized limit of the mean field theory. Both 

length and time scales are shown to diverge at the transition point and 
it has been checked that the nature of the time scale divergence agrees 
Cd ' well with the result obtained from the numerical solution of the mean 

field equation of motion. Similar growths in length and time scales 
r^ . are also observed, as one approaches the transition point, using Monte 

Ch ! Carlo simulations. However, these are not of the same nature as the 

O ■ mean field case. Nucleation theory provides a qualitative argument 

which explains the nature of the time scale growth. To study the na- 
ture of growth of the characteristic length scale, we have looked at 
/\ ' the cluster size distribution of the reversed spin domains and defined 

j^ ■ a pseudo-correlation length which has been observed to grow at the 

phase boundary of the transition. 



The dynamic response of pure Ising systems to time dependent magnetic 
fields is currently being studied intensively JIl. In particular, the response 
of Ising systems to pulsed fields has recently been investigated [0, |^, ^ . The 
pulse can be either "positive" or "negative". At temperatures T below the 
critical temperature T^ of the corresponding static case (without any external 
field), the majority of the spins orient themselves along a particular direction 
giving rise to the prevalent order. In the following, we denote by positive 
(or negative) pulse an external field pulse applied along (or opposite) the 
direction of the existing order. The effects of a positive pulse can be analyzed 
by extending appropriately the finite size scaling technique to this finite time 
window case , and it does not involve any new transition or introduce any 
new thermodynamic scale into the problem. The negative field pulse, on 
the other hand, induces a new dynamic "magnetization-reversal" transition, 
involving completely new length and time scales 0, 0. In fact, we believe, 
the spontaneously occurring dynamic symmetry-breaking transition in Ising 
models under (high frequency) external oscillating fields [^ H occurs actually 
during this "negative" pulse period (and not during the "positive" pulse 
period; compared to the instantaneous existing order in the system) , and the 
universality classes of these two transitions are identical. We report here the 
results of an investigation on the nature of the characteristic length and time 
scales involved in this dynamic magnetization-reversal transition in an Ising 
model under the negative pulsed field. 

In the absence of any symmetry breaking field, for temperatures below 
the critical temperature of the corresponding static case (T < Tc) , there are 
two equivalent free energy minima with average magnetizations +mQ and 
—mo- If in the ordered state the equilibrium magnetization is +mo (say) and 
a very weak pulse is apphed in the direction opposite to the existing order, 
then temporarily during the pulse period the free energy minimum with 
magnetization —rriQ will be brought down compared to that with +mo. If 
this asymmetry is made permanent, then any non- vanishing field (strength), 
which is responsible for the asymmetry, would eventually induce a transition 
from +mo to —mo (in the limit of vanishing field strength). Instead, if the 
field is applied in the form of a pulse, the asymmetry in the free energy wells 
is removed after a finite period of time. In that case, the point of interest 
lies in the combination of the pulse height or strength (hp) and its width or 
duration (At) that can give rise to the transition from +mQ to —mo. We call 
this a magnetization-reversal transition. A crucial point about the transition 




is that it is not necessary that the system should attain its final equilibrium 
magnetization —mo during the presence of the pulse; the combination of hp 
and At should be such that the final equilibrium state is attained at any 
subsequent time, even a long time after the pulse is withdrawn (see Fig. 
1) . The "phase boundary", giving the minimal combination of hp and At 
necessary for the transition, depends on the temperature. As T — ^ T^, the 
magnetization reversal transition occurs at lower values of hp and/or At and 
the transition disappears at T > Tc. 

In the present paper we present an argument that this dynamic transition 
corresponds to infinite time and length scales, all along the phase boundary 
in the hp — At plane at any temperature at T < Tc. We show that the 
relaxation time r and the correlation length ^ both diverge as one approaches 
the phase boundary. In the mean field case, we show (using equations of 
motion linearized in the magnetization) that 

r ~ In ( ) and ^ ^ 

where m^ is the "order-parameter" for the transition, given by the magnetiza- 
tion at the time of withdrawal of the pulse, starting from mo, the equilibrium 
magnetization at the temperature T(< Tc) (see Fig. 1). It may be noted 
that rriwiT, hp, At) = at the phase boundary of the magnetization- reversal 
transition. 

We also show that ^ and r grow sharply as one approaches the phase 
boundary in the Monte Carlo case as well, although the nature of the growths 
are different from the mean field case. We also study the shapes and sizes 
of the reversed spin domains as one approaches the spin-reversal transition 
phase boundary in the Monte Carlo case. We compare the observed growth 
in the relaxation time in this case with that predicted by the nucleation 
theory. 

The Ising model in the presence of an external magnetic field is described 
by the Hamiltonian 

H = --J2 JijSiSj - Y^ hiSi, (1) 

(ij) i 

where Si denotes the spin at ith site, Jij is the cooperative interaction be- 
tween the spins at sites i and j and (...) denotes the nearest-neighbour pairs. 



Here hi is the external field, allowed to be time dependent, and also site- 
dependent to allow investigation of separation-dependent correlations. The 
free energy of the system in the Bragg- Williams approximation is given by[B| 



^ Jijrriimi 



2 Z^ ■JijiiHiii'] 
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(2) 

with TTii = (Si), where (...) denotes the thermal average. In the presence of 
a time and site-dependent field, the time dependent magnetization satisfies 
the Langevin equation 



drrii 
"df 



X 6F 
T 6mi 



X 



Y,Kijmj{t) + 



hi{T) l^^/ l + m,(t) ^ 
T 2 ""yi-m^it)^ 



(3) 



where Kij = Jij/T and A is a constant. Differentiation with the space and 
time dependent magnetic field hi{t) generates the space and time dependent 
susceptibility. After the differentiation we can set hi{t) = h{t) to obtain 
results for a pulsed field uniform in space. Then mj(t) -^ m{t) gives 



dm{t) 
dt 



X 



K{0)m{t) + 



h{t) l^^/ l+m(t) ^ 
T 2 ^\l-m{t)^ 



The resulting equation for the susceptibility, in the Fourier space, is 



dXq{i) 
dt 



X 



K{q)- 



X,it) + ^5{t-t'). 



(4) 



(5) 



l-m?{t) 

Here, K{q) is the Fourier transform of Kij] for small g, K{q) ~ i^(0)(l — g^); 
in the mean field theory -ft'(O) = T^/T. Using @) and (H), we can write 



(6) 



dx,{t) 


'K{q) - 


l-m2(i)J Xq{t) 


dm{t) 


K{0)m{t) - 


hp 1 1 ri+m(i)l 



In the limit when m{t) is small, retaining up to the linear term in m{t), 

dxgjt) _ [K{q) - l]x,{t) 
dm{t) [K{Q)-l]m{t)-Y 
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(7) 



This equation can now be solved in the three different time zones (Fig. 1): 
namely, in the equilibrium regime before the application of the pulse where 
m = ttiq (regime I), the (nonequilibrium) pulsed period regime, at the end 
of which m = 111^ (regime II), and the regime after the pulse is withdrawn 
(regime III) when the system eventually returns to equilibrium (with m{t ^ 
cxd) = —mo if the transition occurs, or = rriQ if it does not). Hence in regime 
II and III, we get the non-equilibrium susceptibility Xq as a function of m(t). 
The solution of (^) also gives the non-equilibrium magnetization m{t), and 
hence we can also arrive at Xg(t). Noting that Xq{t) = Xg when m{t) = mo, 
at the start of regime II, where mo and Xq ^''^^ equilibrium values of the 
magnetization and susceptibility respectively, we can integrate (|^ in that 
regime to obtain 



Xq{t) 



A-a 



where 
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Also integrating the linearized version of (HI) in region II, one gets 



(8) 
(8a) 

m 



m(t) =r + (mo- r)exp[A6(t- to)], (9) 

where b = K[0) — 1. At the end of region II, the value of magnetization is 
given by 



The eqn. 



m^, = m(to + At) = r -|- (mo — F) exp(A6At). 
can therefore be written as 



Xq{t) 



Ac 



exp{Xbaqt) = exp[{K{q) — l)Xt]. 



(10) 



(11) 



In regime III, however, h{t) = and the (initial) boundary condition is 
m(to + At) = m^. Integrating (|^) in this regime one gets 



Xqjt) 
Xq{t0 + At) 



m 



m,. 



(t) 



or 



X,(t)=X^exp[A(ir(g)-l)(t + At)] 



m{t) 



m„ 



(12) 



where use has been made of the eqn. (^. Concentrating on the dominating 
g-dependence of the susceptibihty, one can write 



2f2i 



where the correlation length ^ is defined as 

ln(l/m^ 



^ = ^{m^ 



1 - T/T, 



(13) 



(14) 



This is one of the principal results of this paper, and it shows that the 
characteristic length ^ diverges as the order parameter m^ goes to zero. 

Consider now the t dependence arising in Xq=o(t) through the factor 
m{tYt. Solving (^) in regime III yields 



m{t) = rriyj exp[A6{t — (to + A^)}], 



(15) 



which shows that long time is required to attain moderate values of m{t) 
starting from low values of m^. Especially, starting from time t = to + ^'t, 
the time taken by the system to reach the final equilibrium value is defined 
as the relaxation time r of the system. Therefore from (|T5]) we can write 



r 



T 



X\Z 



/mo 
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(16) 



The growth of the time scale occurs in Xg=o(^) too through the m(t) depen- 
dence : 

^^ ~exp[A6{t-(to + At)}]. 



Xq=0it) 



m„ 



Eqn. (|14D and (|16]) can be used to establish a relationship between r and ^ : 



In ( — ) ~ ^e- 



(17) 



This corresponds to critical slowing, with the characteristic time diverging 
with the characteristic length with the dynamical critical exponent 

z = 2. 

The above results are obtained in the linearized limit of the mean field 
eqns. of motion @) and (H). We also measured, solving the full dynamical 
equation (^ numerically, the relaxation time r by computing the time re- 
quired by m{t) to reach the final equilibrium value ±mo, with an accuracy 
of O(10~^), from the time of withdrawal of the pulse (in regime III). Fig. 
2 shows that this r indeed diverges as one approaches the phase boundary, 
where m^ = 0. In fact, the numerical results are observed to fit very well 
with the analytic result ( p^ (shown by the solid line in Fig. 2). 

The divergence of both the time and length scale were also investigated 
at low temperatures by employing Monte Carlo methods. Simulations on 
a square lattice of typical size L = 200 with periodic boundary conditions 
indicated an exponential growth of the time scale : 

T ~ exp[-c(T) \m^\], (18) 

where c{T) is a constant depending on temperature only. Further, finite size 
scaling of the order parameter relationship 

m^~ \K-h;f (19) 

is consistently with (3 = 0.90 ±0.02 and with a correlation length divergence 
with z/ = 1.5 ± 0.3. (Here h^ is the critical value of the pulse field hp, 
making m^ = at the end of regime II). These results qualitatively compare 
with the divergence of scales at the transition point predicted by the mean 
field treatment. However, the growths of the time and length scales are 
quantitatively of different nature to that of the mean field case, because at 
low temperatures droplet growth is a dominant mechanism. The growth of 
droplets of size / is associated with an activation energy^] E{1) = —2hpl'^ + 
al'^~^, where a is the surface tension. Using the relationship between / and hp 
at the energy minimum together with (18) at small rriw gives a characteristic 
time 

exp[-ci(r) I m„ I'/" (A=)'-=]. (20) 



f'i- 



r ~ exp 
Since (3 is close to unity, this is consistent with the observed relation 



The typical size of cluster or domain of reversed spins provides a qual- 
itative idea about the correlation length of the system. In order to study 
the growth of the typical reversed- spin domain size, we define a pseudo- 
correlation length C, as follows : 

where Ug is the number of domains or clusters of size s and the radius of 
gyration Rg is defined as R^ = J2i=i I ^i — vq p /s, where r^ is the position 
vector of the ith spin of the cluster and tq = J2i=i{ri/s) is defined as the 
centre of mass of the particular cluster. The pseudo-correlation length ^ is 
observed to grow to system size order as one approaches the phase boundary 
(Fig. 3); thereby providing further indication of the growth of a length scale. 
It should be noted that, as in the static transition in the pure Ising system, 
the length ^ is distinct from the correlation length . 

In the linear limit of the mean field dynamics, it has been possible to 
show the divergence of both the length and time scales at the magnetization- 
reversal transition phase boundary. Sharp growth of these scales has also 
been observed in the Monte Carlo case, studied in two dimension. Here, 
we looked at the size distribution of the clusters or domains of reversed 
spins whose average size was observed to grow at the phase boundary of the 
transition. 
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Figure Captions 

FIG. 1. Schematic time variation of the pulsed field h{t) and the corre- 
sponding response magnetization m{t) for two different cases. The solid line 
indicates no magnetization-reversal case whereas the dashed line indicates a 
magnetization- reversal. 



FIG. 2. Divergence of the relaxation time in mean field limit for T = 0.8 and 
At = 20 (from numerical solution of eqn. (4)). The solid line indicates the 
corresponding analytical estimate (eqn. (15)). 



FIG. 3. Growth of the pseudo-correlation length ^ for different system sizes 
in the Monte Carlo study on a square lattice of size L x L. 
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